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We study the ground states of one-dimensional Bose-Bose mixtures under harmonic confinement. As we vary 
the inter-species coupling strength up to the limit of infinite repulsion, we observe a generalized, composite- 
fermionization crossover. The initially coexisting phases demix as a whole (for weak intra-species interactions) 
and separate on an atomic level (for strong intra-species repulsion). By symmetry, the two components end 
up with strongly overlapping profiles, albeit sensitive to symmetry-breaking perturbations. Different pathways 
emerge in case the two components have different atom numbers, different intra-species interactions, or different 
masses and/or trap frequencies. 
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I. INTRODUCTION 

The availability of cold atoms has made it possible to real- 
ize many fundamental quantum systems. Building on the sem- 
inal realization of Bose-Einstein condensation (H2l, mixtures 
composed of, say, two different atomic species have come into 
the research focus. Aside from Bose-Fermi yl |4J or Fermi- 
Fermi mixtures [5], whose potential for studying phenomena 
as diverse as impurity effects or superconductivity has been 
recognized more recently, two-component bosonic mixtures 
have received much experimental 
oretical attention (see 



10|] and the- 
m\M and 

Refs. therein). The interplay between intra- and inter-species 
forces gives rise to many effects not accessible with single- 
component Bose gases, including phase separation and mod- 
ified superfluid-insulator transitions IU 4) , M 5I1. quantum emul- 
sions 1 16], and spin-charge separation HISTI . 

Most studies so far have focused on the regime of rela- 
tively weak interactions, where the physics can be described 
well in terms of mean-field or — in lattice geometries — simple 
lowest-band models. However, interatomic forces can be ex- 
perimentally tuned to a large extent via Feshbach resonances 
11911 . In particular, in quasi-one-dimensional systems, which 
emerge under strong transversal confinement, it is possible to 
exploit confinement-induced resonances B20I1 to explore the 
regime of strong correlations I2H l22tl . For infinitely repulsive 
bosons, this is known as the fermionization limit, in allusion 
to the fact that the system can be mapped exactly to an ideal 
Fermi gas [23fl. Here the exclusion principle in a sense emu- 
lates the effect of hard-core interactions, to the extent that the 
bosons share local aspects with their fermionic counterparts, 
whereas nonlocal properties such as their coherence and mo- 
mentum distribution are very different. The basic crossover 
from the weakly interacting trapped Bose gas to the fermion- 
ization limit had been predicted from a thermodynamic-limit 
perspective |[24l [25tl and interpreted in terms of a mean-field 
picture Ir26ll . By contrast, it is only recently that its micro- 
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scopic mechanism has been investigated within an ab-initio 
framework JH, M, S HI • 

In this work, we tackle the obvious question of how the 
fermionization crossover for the one-component Bose gas ex- 
tends to a trapped two-component mixture. By way of anal- 
ogy, tuning the inter-component coupling strength to the in- 
finitely repulsive regime (for fixed intra-species interactions) 
may be regarded as composite fermionization. Here, a re- 
cent study has extended the standard fermionization map to 
mixtures of two identical particle species with both intra- and 
inter-species hard-core interactions [31]. Apart from this spe- 
cial borderline case, little is known except for a classification 
of the low-energy modes in the harmonic-fluid approximation 
ifTTll . Here we study the crossover from weak to strongly re- 
pulsive couplings between two components under harmonic 
confinement. We will show that this composite fermioniza- 
tion can lead to demixing, and lay out how it depends on the 
intra-species interactions, on the densities of the two compo- 
nents, as well as on the masses and trapping parameters of 
each species. 

Our paper is organized as follows. Section [TT] introduces 
the model and briefly reviews the fermionization map and its 
extension to mixtures. In Sec. [HI] we give a concise presen- 
tation of the computational method. Section [TV] first explores 
the completely symmetric setup, where both components have 
equal atom numbers, interaction constants, masses, and see 
the same harmonic trap. The subsequent Sec.[IV]in turn shows 
what different phase-separation scenarios emerge if these con- 
straints are relaxed one by one. 



II. THEORETICAL BACKGROUND 
A. Model 

We consider a mixture of two distinguishable bosonic 
species, which we shall label "A" and "B". These may corre- 
spond to atoms with unequal nucleon numbers — be it different 
isotopes or altogether different species — or possibly different 
hyperfine states of one and the same species. Furthermore, 
we assume these to be confined to quasi-one dimension (ID), 
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such that the transverse degrees of freedom may be integrated 
out. The effective low-energy Hamiltonian for an arbitrary 
mixture of N — N A + Nb atoms then reads 



H= h *+ h ab, 



<t=A,B 



where the single-species Hamiltonian H a and the inter- 
species coupling Hab read 



H„ 



E 

i=l 



2M„ 



N A N b 

Hab = ^ ^ ffAB^(^A,q - X B .b)- 
a=l b=l 



Here we consider harmonic trapping potentials U a (x) = 

tor units a a 
1 by exploiting 



. By rescaling to harmonic-oscillator units <za 



\J h/MpjaJA, one can eliminate Ma 
the scaling 



UJA 



H A (M A , w A , g A ; X a ) = fHj A H A (l, 1, g' A ; X' A ), 



with X' A = 



( x A,i' 



' X A,N A 



X A /a A and # A = 



g A ^/M A /h 3 uj A . 

For numerical reasons, we regularize the delta-function in- 
teraction by a normalized Gaussian of width much smaller 
than the inter-particle distance; see Ref. [28] for details. 



B. Fermionization 

The (single-component) ID Bose gas has the peculiar prop- 
erty that it is isomorphic to a system of identical fermions. 
In particular, the standard Bose-Fermi map relates the many- 
body wave function of hard-core bosons (obeying the bound- 
ary condition \?| . = 0, i < j, which corresponds to 
taking the ID interaction strength g — * oo) to that of non- 
interacting fermions \?_: 



A(X) := J|sgn(a;j - Xj). 

i<j 



Specifically, the ground state is given simply by the absolute 
value of the non-interacting fermionic ground state, vJ/ = 

| \3> . o I - This makes it tempting to think of Pauli's exclusion 

principle as emulating the effect of the repulsive interactions 
(or vice versa), which is why the limit g — > oo is commonly 
referred to as fermionization. Note that, since A 2 — 1, all 
local quantities will coincide with that computed from the 
fermion state. Specifically, this is the case for the density 
Pn = 1^1 2 and any derived quantities, such as the reduced 
(one- or two-body) densities. However, nonlocal quantities 
such as the momentum distribution, may differ dramatically 
from the fermionic ones. 

The standard Bose-Fermi map above has recently been ex- 
tended to mixtures of two identical species A = B (i.e., equal 
masses and potentials) with hard-core intra- and inter-species 



interactions, g a = gAB — ► oo. Its wave function is trans- 
formed to that of a system of N — J2 a identical fermions 
lITill . For the special case of A and B being bosonic, this gen- 
eralized Bose-Fermi map \P = A^> _ reads 

A(X A ,X B ) = Aa(Xa)Ab(Xb)Aab(Xa,Xb), (1) 



where A a (X a ) = Jli 



<i<j<N c 



sgn(x CT) , 



is the stan- 



dard map restricted to subsystem a, and Aab(^a,-^b) 



Il^i Ylb=i sgn(a;A,a — %~B,b) serves to impose hard-core 
boundary conditions on inter-species collision points. In the 
case of harmonic trapping, where the single-particle orbitals 
are known analytically, the solution may even be written down 
explicitly HI 

(X A ,X B )) oce-l*' 2 / 2 [] \ x i- x A- ( 2 ) 

l<i<j<N 

IIL COMPUTATIONAL METHOD 

Our approach relies on the numerically exact Multi- 
Configuration Time-Dependent Hartree method J32I l33ll . a 
quantum-dynamics tool which has been applied successfully 
to systems of few identical bosons (see ll28l I29I [34l l35l l36ll ). 
Its principal idea is to solve the time-dependent Schrodinger 
equation = H^(t) as an initial-value problem by ex- 

panding the solution in terms of direct (or Hartree) products 



rN B 



(g) if 



(N). 

3N 



*(o=x;^(*) $ j(*)- 



(3) 



The (unknown) single-particle functions ip^ (j = 1, . . . , n K ) 
are in turn represented in a fixed primitive basis implemented 
on a grid. In our case of, where particles of each species are in- 
distinguishable, the single-particle functions within each sub- 
set k € {1, . . . , N A } and {N A + 1, . . . , N} are of course 

identical (i.e., we have {cp^}, with j a < n a ). This, along 
with the correct symmetrization of the expansion coefficients 
Aj, ensures permutation symmetry within each subset A,B. 

Note that in the above expansion not only the coefficients 
Aj but also the single particle functions ipj are time depen- 
dent. Using the Dirac-Frenkel variational principle, one can 
derive equations of motion for both Aj, tpj i33ll . Integrating 
this differential-equation system allows us to obtain the time 
evolution of the system via (01. This has the advantage that the 
basis {$,/(£)} is variationally optimal at each time t. Thus it 
can be kept relatively small, rendering the procedure very ef- 
ficient. 

Although designed for time-dependent simulations, it is 
also possible to apply this approach to stationary states. This 
is done via the so-called relaxation method IB7I1 . The key idea 
is to propagate some wave function ^(0) by the non-unitary 
e~ Hr (propagation in imaginary time.) As r — ► 00, this expo- 
nentially damps out any contribution but that stemming from 
the true ground state like e~( Em ~ E °> T . In practice, one relies 
on a more sophisticated scheme termed improved relaxation 
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l38tl . which is much more robust especially for excitations. 
Here (^\H is minimized with respect to both the coeffi- 
cients A j and the orbitals (pj . The effective eigenvalue prob- 
lems thus obtained are then solved iteratively by first solving 
for Aj with fixed orbitals and then 'optimizing' ipj by propa- 
gating them in imaginary time over a short period. That cycle 
will then be repeated. 

IV. COMPOSITE-FERMIONIZATION TRANSITION 

In contrast to the case of a single bosonic species, binary 
mixtures offer a plethora of different parameters, making the 
physics richer and less straightforward: In principle, we may 
have different atom numbers N = Na + Nq, different masses 
M ct= a,b. intra- and inter-species couplings g a (<?abX an d 
species-dependent traps U a (x). In this section, in order to 
illustrate the basic mechanism of the crossover from weak to 
strongly repulsive inter-species interactions, <7ab € [0,oo), 
we focus on the simplest, symmetric setup where 

N a = —■ M a = 1; g a = g; U a (x) = ^x 2 (a = A, B). 

In this case, H has an exact permutation symmetry between 
species A and B. This idealized situation may correspond to 
two internal states of the same species or, ignoring slight mass 
deviations, two different isotopes, where <?ab is tuned via the 
inter-species scattering length. Actually, in the special case 
where g a — qab = g, this system maps to a one-component 
Bose gas with N = Yl a N a atoms ll28l I29H30I1 (for any num- 
ber of components and any N a , for that matter) - up to permu- 
tational degeneracies, which are not that severe for the ground 
state. 

Here, by contrast, we are interest in the following ques- 
tion: What happens to the mixture — g a being fixed — when 
<7ab is varied up to the hard-core limit? (This we refer to 
as composite fermionization of the subsystems A and B, de- 
spite the general lack of a Bose-Fermi mapping as Eq.Q]) To 
get an impression of that crossover, let us start with the case 
of two almost ideal Bose gases, g a = 0.4, each consisting 
of N a = 2 atoms (similar results hold for larger atom num- 
bers). Figure Q] displays the evolution of the density profile 
p(x) = p a (x), measuring the probability distribution for find- 
ing one a atom at position x. Obviously, for $ab — > 0, the 
total state = \I> A ® \]> B simply consists of two uncorrelated 
"condensates" (ty a = 4>f " for g = 0), slightly smeared out 
due to repulsion. Increasing gAB leads to an ever deeper dip 
in the profiles. This should be contrasted with the case of two 
single fermionized bosons, N a — 1 [39]- The dip in Fig.[T]is 
much more pronounced, which is indicative of phase separa- 
tion, if symmetry screened: pa — Pb are completely identical 
by symmetry. However, this only corresponds to an ensem- 
ble average - in a single measurement, we will always find all 
Na atoms on one side of the trap and Nb atoms on the other. 
This claim is underscored by Fig. \V[bottom), which reveals 
the evolution of the two-body densities p a .a'- If we were to 
measure, say, the first A-type boson at xa,i ~ 1, then we are 
sure to find the second A boson also in that region xa : 2 ~ 1 
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Figure 1: Composite fermionization of a mixture with iVo-— a.b = 2 
bosons with intra-component interaction g„ — 0.4. Top: density 
profiles p(x); Bottom: Two-body correlation functions p a <j (xi, xz) 
and pab(xa, xb) for inter-species couplings (?ab = 0.4, 4.7 and 25 
(from left to right). 



and not on the left, and vice versa. By contrast, the probabil- 
ity for subsequently finding a B particle at the same position 
is virtually zero, as dictated by the hard-core boundary condi- 
tion, ^I^a a =x B b = (ffAB — ► oo). This makes it tempting to 
think of this an entangled state of the form 

\N A , 0) ® |0, N B ) + |0, N A ) ® |JV B , 0), 

where |nL,TiR,) S H CT denotes a state with njj,R atoms lo- 
calized on the left (right). It should be noted that, even for 
g a = 0, there is no simple mapping to fermions as (fl}, since 
the hard-core condition is imposed only on inter-species colli- 
sion points, and thus the information about which fragment 
the individual coordinates belong to needs to be retained. 
However, in our special case of a harmonic trap, it is natural 
to conjecture that the exact solution is given by a modification 
of the single-species fermionization limit 114011 . 

* gAB ^oo(X) OC e - |X|2/2 Y[ I^A.o ~XB,b\, 
a<N A ,b<N B 

which obeys the correct boundary conditions at points of inter- 
species collisions. Trusting that logic, an analogous extension 
should hold for the homogeneous system [23]. 
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Figure 2: Same as Fig.Q] but with qa = <?b = 4.7. 



A similar pathway is encountered for two more strongly 
interacting components, g a — 4.7 (see Fig. [2] top). At 
9ab = 0.4, we have more or less two uncorrelated clouds, 
which are governed by the desire to reduce their infra-species 
interaction energy. As <?ab = 4.7 reaches g a , this turns into 
a trade-off between avoiding the own species just as much as 
the other component. Letting <?ab — > oo, the inter-particle 
repulsion takes over, and a similar phase-separation tendency 
of A and B as before may be recognized in Fig. [2] In contrast 
to the "condensate" case, however, the separation of the two 
peaks is not pronounced as each hump is quite smeared out 
in itself due to the intra-species repulsion. This is illuminated 
further by the two-body densities (Pig.[2j: Here the pattern for 
p aa (xi ) X2) at <?ab = 25 is modulated by a correlation hole 
at x\ = X2 due to intra-species repulsion, as compared to the 
weakly interacting components (Fig.[T|). This explains the two 
broadened peaks in p a (x). 

So far, we have seen that the components tend to separate 
when the inter-species repulsion overwhelms the intra-species 
one. This naturally brings up the question of the fate of two 
initially fermionized components, as shown for g a = 25 in 
Fig. [3] Notably, by the conventional Bose-Fermi map, this re- 
lates to a Fermi-Fermi mixture. Weak couplings <?ab = 0.4 
pass the two fermionized clouds largely unnoticed, which ex- 
hibit N a — 2 characteristic humps in p a {x) [40]. However, 
for larger values $ab =4.7, the profiles slowly rearrange to 
a more complex structure, which culminates in a profile with 
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Figure 3: Same as Fig.Q] but with gA = <7b = 25. 



N = 4 wiggles at (?ab = 9a = 25. The density oscillations 
signify that each of the four atoms seeks an isolated spot, ir- 
respective of its species. That interpretation is supported by 
the plots of the two-body densities p aa = pab, which for 
<7ab = 25 reveal the checkerboard pattern familiar from the 
single-boson crossover ll28ll . This should be contrasted with 
the intermediate regime where ^ab = 4.7 < g a : Here two, 
say, A atoms are still localized on the left and on the right side 
as for 5ab = 0. Upon measuring an A atom at, say, xa ~ 1.5, 
the two B atoms will likely be found at either xb ~ 0.5 or 
ib ~ —1.5, this way remaining isolated from each other but 
also avoiding the A atom. 

Note that, in agreement with our earlier remarks, the case 
9a = 9ab relates to a single-component Bose gas, which in 
turn maps to an ideal Fermi gas via (Q~|i in the limit #ab — * oo. 
As in that case, for N 3> 1 these N peaks become ever tinier 
modulations on the envelope density, which for a harmonic 
trap can be computed as p(x) = y/2N — x 2 /Ntt l4lll . 

At this stage, we should point out that this limit is highly de- 
generate: For one thing, there is a permutation degeneracy be- 
tween A and B particles. Second, in the limit <?ab — * oo, the 
ground-state wave function (which is non-negative) degener- 
ates with the fermionic one by the Bose-Fermi map and, since 
no specific permutation symmetry is imposed when treating 
the two components as distinguishable, all solutions even with 
mixed A-B-exchange statistics are permissible 113 ill - 
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Figure 4: (color online) Instability under symmetry-breaking pertur- 
bation, U a (x) = \x 2 + dcrX (cLa = — g(b = 0.1). Left: Densities 
pcr(x) for g a = 0.4; Right: g a — 25. Shown are the coupling 
strengths jab = 0.4, 4.7, 25 from top to bottom. 



Symmetry-breaking instability. By permutation symme- 
try of H between A and B, the density profiles p a are iden- 
tical, even in the limit gAB 3> ga, and thus trivially cannot 
exhibit phase separation. However, if we deal with two dif- 
ferent species, it is conceivable that these two feel slightly 
different trapping potentials, where the deviations are much 
weaker than the mean trapping and inter-particle forces but 
only serve to break the symmetry. In particular, imagine that 
U a (x) — ix 2 + d a x, such that the trap centers be shifted 
by c£a = — ds -C 1, see Fig. [4] Expectedly, for weak cou- 
plings gAB, the profiles are barely affected. However, toward 
stronger inter-species repulsion, this tiny perturbation is the 
last straw needed to make the two phases demix completely. 
Similar results hold also for different densities, Na ^ Nb- 
This makes it even more inviting to think of the symmetric 
profiles in Figs. Q~|-[3]as averages over the equivalent configu- 
rations with A (B) being on the left (right), and the other way 
around. 



V. PHASE-SEPARATION SCENARIOS 

So far, we have studied a completely symmetric setup, 
where only the inter-species interactions were permitted to 
differ. This way, a symmetry-breaking perturbation was 
needed to reveal the hidden phase separation. Although not 
experimentally unrealistic, this scenario is somewhat artifi- 
cial. We now want to relax the above symmetry constraints 
step by step and discuss the wealth of different demixing path- 



ways if the two components have different particle numbers 
(Sec. IV At , different internal interaction strengths (Sec. IV Bt , 
and have different masses and/or trap frequencies (Sec. IV CI . 



A. Density-assisted demixing 

An obvious question regarding our findings in the previous 
section is: What happens in the case of unequal particle num- 
bers, A*a 7^ -Mb? Figure [5] illustrates this on the example of 
two weakly interacting components, g a = 0.4 (left column), 
where Na = 3 is larger than Nb = 2. As the Bose-Bose 
coupling gets stronger, </ab = 1-3, one observes that the low- 
density phase B moves to the outer edge, thus "sandwiching" 
high-density A component in the middle. This well-known 
phenomenon traces back to the fact that the coupling en- 
ergy (Fab) = N A N B gAB J pab{x,x)(1x scales with N A N B 
compared to the individual energies (H a ) oc N a . Thus for 
the smaller B component, it is less expensive to move to the 
higher potential regions. Not unexpectedly, we have found 
this to be even more pronounced for Na S> Nb- Note that, 
close to the composite-fermionization limit ((?ab = 4.7, 25), 
both components develop two humps in the density, if much 
more pronounced for B. This is indicative of a superposition 
state similar to that in Sec. [IVj The B atoms are found on the 
right and the A atoms on the left, and vice versa, only that the 
shift for A is much smaller due to their higher density. 

A similar pathway again exists for non-negligible intra- 
component interactions, g a — 4.7 (Fig. |5}. As qab — ► oo, 
the initially mixed phases separate: The profile pa(x) devel- 
ops a clear-cut peak at x = 0, whereas B is again driven to 
the boundary. Even though, on the face of it, this looks dif- 
ferent from the weakly interacting case, this density pattern 
can be understood in complete analogy: The two components 
are isolated on the left and on the right, respectively; however, 
due to the larger atom number Na and the repulsion pressure, 
A tends to be more in the center on average. 

An entirely different situation is encountered in the "Fermi- 
Fermi"-like setup with g a = 25 (Fig. [5]). For intermediate 
gAB = 15 < g a , a phase forms where the A atoms local- 
ize at three discrete spots such that the two B atoms fill the 
two holes in between. For <?ab = 25, by contrast, the N 
atoms completely localize atom by atom just like in the case 
of a 2+2 mixture, in agreement with the extended Bose-Fermi 
map (fl}. In analogy to Sec. lIVI these will demix under slight 
symmetry-breaking perturbations into one phase with Na = 3 
density wiggles on, say, the left side, and Nb = 2 on the right. 



B. Interaction-assisted demixing 

Up until now, we have assumed comparable interactions 
within each component. Of course, it is of fundamental con- 
cern what the composite-fermionization crossover looks like 
in the case where one species is more strongly repulsive, in- 
cluding as a special case a "Bose-Fermi"-type mixture of one 
weakly interacting and another, fermionized component. 
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Figure 5: Demixing for different particle numbers Na = 3, Nb = 2: Density profiles pa(x) (top) and Pb(x) (bottom) for intra-species 
interaction strengths g a = 0.4, 4.7 and 25 (from left to right). 
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1. For <7ab < 9A, the weakly interacting central B cloud 
is barely affected; at the same time, the strongly inter- 
acting A bosons move slightly to the outside, thus cut- 
ting down on both intra- and inter-species interaction 
energy. 

2. By contrast, for g&B = 25 > gA, this partial separation 
is no longer enough: Now the B cloud splits up, sig- 
nifying the formation of the entangled state discussed 
before with Na atoms on the left and TVb atoms on the 
right, and vice versa. Note that, owing to the strong 
inter-species repulsion in A, the two humps in pa are 
washed out strongly. 

We stress that only regime (1.) exists for a Bose-Fermi-type 
mixture, i.e., where qa — * oo: The minimum-energy state for 
infinitely large (?ab then has all B atoms in the center and A 
on the edges. 

Coherence aspects. At this point, it is worthwhile 
dwelling for a moment on the coherence properties of bosonic 
mixtures, as reflected in the reduced one-body density matrix 



p a {x,x') = 

distribution 



(x\pa 



Ml 



and, closely related, the momentum 



Figure 6: Demixing for different intra-component interactions gA = 
4.7, gs = 0.4: Density profiles pa(x) (top) and pb{x) (bottom) for 
an Nc = 2 mixture at different couplings #ab = 0.4, 4.7 and 25. 



An illustrative example is given in Fig. [6] displaying the 
composite-fermionization crossover for an N a — 2 mixture 
with gA = 4.7 > gB — 0.4. We distinguish two regimes: 



p a (k) = 2Tr(k\p£ ) \k) = dx dx'e- tk{x - x ">p a (x,x') 



It has been demonstrated for identical bosons [29, 30] how, 
in the course of fermionization, the zero-momentum peak 
p(k = 0) — related to the fraction of condensed bosons — is at- 
tenuated and redistributed toward higher momenta, culminat- 
ing in a characteristic decay p(k) k ~^° ck~ A as predicted for 
hard-core short-range interactions 14211 . Equivalently, the off- 
diagonal long-range order, measured by p±(x, —x) as x — > oo, 
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Figure 7: Coherence properties of the system in Fig. [6] Top: mo- 
mentum distributions p a (k); Bottom: One-body density matrices 
Pa(x a , x'a) for gAB = 0.4, 4.7 and 25 (from left to right). 



is strongly reduced. We generally find the same two mech- 
anisms at work here, which we exemplify in Fig. [7] For 
stronger inter-species repulsion, the high-momentum tail in 
p a (k) becomes more pronounced. Interestingly, for the com- 
ponent B with weaker interaction, the k = peak starts di- 
minishing right away, while the strongly repulsive A compo- 
nent first sees a sharpening at zero momentum (<?ab = 4.7). 
This derives from its initial derealization so as to move away 
from B [cf. pa(xa, x' a ) in Fig. [TJ^ which allows the A atoms 
to spend less kinetic energy (hp 2 ) oc (Afc) 2 . 



C. Trap-induced demixing 

After having explored the effect of different densities or 
interaction strengths on the composite-fermionization path- 
way, let us now relax the condition of equal masses and trap- 
ping potentials (here: frequencies). In this case, the sys- 
tem no longer maps to a single-component Bose gas even for 

9 a = 3AB- 



1. Different confinement lengths 

Assume that we have a nontrivial mass ratio, i.e., Mb > 1 
without loss of generality, with an otherwise symmetric pa- 
rameter set. The effective oscillator length of the B atoms will 
then be reduced by a factor of ob = 1/ a/Mb < 1- This situ- 
ation is visualized in Fig. [8] for the choice Mb = 9. At weak 
couplings, Pb( x ) is simply constricted at the trap center, while 
Pa(x) extends over a much larger region. As we switch on 
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Figure 8: Demixing for different masses, Mb /Ma = 9: Density 
profiles Pa(x) (top) and Pb(x) (bottom) for a 2 + 2 mixture with 
g a = 0.4 at different couplings #ab. 



the interaction between the components, the B atoms remain 
unmoved, whereas the A bosons are gradually driven toward 
the outside. This is intuitive: The former component roughly 
feels an average Hamiltonian 



tr A [H A Bpf A 



Hb + gABN A ^PA(xB,b), 
6=1 



and likewise for A. Since the heavy B atoms are effectively 
frozen at the center, where pa{xb) ~ Pa(0) changes slowly, 
they only feel a constant energy shift due to the presence of A 
atoms. By contrast, the latter ones see an effective "potential 
barrier" pb(^a) ~ ^(^a) which varies only in a small region 
about zero. 

That phase-separation mechanism is largely insensitive to 
the intra-species interactions g a : We have confirmed these re- 
sults also for, e.g., two quasi-fermionized components. Also 
note that a similar scale separation persist for the case of dif- 
ferent frequencies but equal masses, i.e., cjb/^a 3> 1. The 
different effective interaction felt by B, g' B = 3b\/Mb/wb, 
and the modified energy scale, lob, do not qualitatively alter 
the picture above. 



2. Different energy scales 

Let us now look into the complementary case where the os- 
cillator lengths be equal, as = 1/ x/Mbujb — 1, but such that 
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Figure 9: Demixing for different energy scales, Mb = 1/i^b = 3: 
Density profiles pa(x) (top) and pb(x) (bottom) for a 2 + 2 mixture 
with qa = 25 and <jb = 0.4, at different couplings jab. 

the energy scale lj-q ^ 1 shall be different. In order words, a 
stronger localization by virtue of a larger mass is compensated 
by a shallower trap for the B species. This option is sounded 
out in Fig.|£l where Mb = 3 = cjg 1 . The two profiles still 
overlap for weak couplings <?ab = 0, 1.3. For sufficiently 
strong inter-species repulsion, though, it apparently becomes 
beneficial for the B atoms to spread out to larger x so as to 
segregate from A, which in turn is compressed on the inside. 
This is particularly striking in the setup captured in Fig. [9] 
Here the A component is squeezed even though it is fermion- 
ized and thus possesses a high internal pressure. The reason 
for that counter-intuitive behavior is simply that the potential- 
energy costs for the B phase are lower by lo-q = 1/3. 



VI. CONCLUSION 

We have studied binary few-boson mixtures in a one- 
dimensional harmonic trap throughout the crossover from 
weak coupling to strong inter-component repulsion. Depend- 
ing on the intra-species interactions, different pathways to a 
new kind of "composite fermionization" have shown up: For 
two weakly interacting Bose gases, the two phases segregate 
as a whole, where the demixing for equal densities is obscured 
by symmetry-induced entanglement fragile to displacement of 
the trap. By contrast, for two strongly repulsive components, 
demixing occurs on the atomic level. 

If one component has a lower density, then it tends to de- 
localize toward the outer edge, while the high-density phase 
is compressed in the center. Furthermore, in case one com- 
ponent is far more repulsive, the crossover exhibits an inter- 
mediate regime where that species forms a shell around the 
central, weakly interacting, component; only for large inter- 
species couplings do they fully segregate. This is accompa- 
nied by an increase (decrease) of the central momentum peak 
for the strongly (weakly) interacting species. Finally, for dif- 
ferent mass or frequency ratios, one component freezes at the 
trap center, such that it acts as an effective potential barrier for 
the more mobile species. 

The small mixtures of strongly repulsive atoms studied here 
should be experimentally accessible. The preparation and de- 
tection techniques required are similar to those already avail- 
able for few bosons of a single species. The interaction forces 
may be tuned independently over a wide range by varying 
the (inter- and intra-species) scattering lengths as well as the 
transverse confinement, which parametrically modifies the ef- 
fective one-dimensional coupling strengths. 
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